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Abstract 

Transport of scalar fields in compressible flow is investigated. The effective 
equations governing the transport at scales large compared to those of the advecting 
flow V are derived by using multi-scale techniques. Ballistic transport generally takes 
place when both the solenoidal and the potential components of v do not vanish, 
despite of the fact that v has zero average value. The calculation of the effective 
ballistic velocity Vb is reduced to the solution of one auxiliary equation. An analytic 
expression for Vb is derived in some special instances, i.e. flows depending on a 
single coordinate, random with short correlation times and slightly compressible 
cellular flow. The effective mean velocity Vi, vanishes for velocity fields which are 
either incompressible or potential and time-independent. For generic compressible 
flow, the most general conditions ensuring the absence of ballistic transport are 
isotropy and/or parity invariance. When vanishes (or in the frame of reference 
comoving with velocity Vb)^ standard diffusive transport takes place. It is known 
that diffusion is always enhanced by incompressible flow. On the contrary, we 
show that diffusion is depleted in the presence of time-independent potential flow. 
Trapping effects due to potential wells are responsible for this depletion. For time- 
dependent potential flow or generic compressible flow, transport rates are enhanced 
or depleted depending on the detailed structure of the velocity field. 



1 Introduction 



One of the most interesting issues in statistical mechanics is related to the effects of a 
given microscopic dynamics on phenomena occurring at much larger scales. A well-known 
example is provided by kinetic theory. When the Knudsen number is small, i.e. the scales 
of interest are much larger than the mean free path, hydrodynamic equations can be 
derived. In these equations only the macroscopic scales are involved. The microscopic 
small-scale degrees of freedom are averaged out and do not appear explicitly. This impor- 
tant point is physically related to the fact that macroscopic and microscopic space-time 
scales are widely separated (scale separation). In this situation it is indeed just the mean 
cumulative result of many almost independent effects which is relevant for the large-scale 
dynamics. This also explains the universality of hydrodynamic equations, whose struc- 
ture is essentially dictated by the symmetries and the conservation laws of the microscopic 
dynamics. The details of the dynamics only affect the numerical value of the transport 
coefficients, e.g. the viscosity in Navier-Stokes equations. A remarkable application of 
this universality has been made in lattice gas methods |^. 

The analogy with kinetic theory is very useful for any system involving scale sepa- 
ration. Let us consider for example the transport of a scalar field B by a velocity field 
V [Q. Such a problem arises, for example, studying the spreading of particles advected 
by a given velocity field and subject to molecular diffusion. At space-time scales much 
larger than those of v, we physically expect the dynamics of O to be governed by an 
effective transport equation. From a technical point of view, exploiting the scale separa- 
tion to derive the effective equations involves a singular perturbation problem (see, e.g., 
Ref. 0]). To deal with it, multi-scale techniques (also known as homogenization) have 
been mostly used 0. A well known application of these techniques is for scalar trans- 
port by incompressible velocity fields (divergence- free). In Ref. it was shown that the 
large-scale dynamics of the scalar in the presence of scale separation is indeed governed 
by an effective equation, which is always diffusive. The calculation of the effective diffu- 
sivity is reduced to the solution of one auxiliary equation. From this equation it follows 
that the effective diffusivity is always larger than the molecular one, i.e. incompressible 
flow enhance diffusion. The auxiliary equation can be solved and the effective diffusivity 
calculated exactly in some special instances, e.g. parallel flow 0, or asymptotically in 
the limit of high Peclet numbers for cellular flow P, H . Numerical techniques are other- 



wise generally needed |T^, |Tl|. Variational principles [T^ for the effective diffusivity have 
also been derived both for static and time-dependent incompressible flow. Note that for 
random flow arbitrarily small wavevectors might be excited. The existence of a range 
of scales where scale separation applies is thus not guaranteed. A necessary condition is 



that the variance of the vector potential be finite [12, 14|. When the previous condition is 



violated, anomalous transport might occur. Hereafter and in the sequel we shall suppose 



2 



that, for scales sufficiently large, scale separation holds and we shall be interested in the 
dynamics at these scales. Another transport phenomenon which has attracted much inter- 
est is momentum transport in Navier-Stokes incompressible flow. This system provides 
a very interesting counterexample to the possible belief that, if scale separation holds, 
the dynamics of the large scales is necessarily diffusive. The so-called anisotropic kinetic 
alpha effect generically appears in the absence of parity-invariance and isotropy. Fur- 
thermore, even when these two conditions are satisfied, the effective viscosity might be 
negative 0, |T^. Large scales are then strongly amplified, nonlinearities become relevant 
and the dynamics at large scales presents some peculiar properties |[T^, |TP|. 

The aim of our paper is to investigate scalar transport in compressible flow. The 
results presented here illustrate the important differences between scalar transport by 
compressible and incompressible flow. First, it is known that transport rates are always 
enhanced by incompressible flow. We show that for generic compressible flow this prop- 
erty is not valid. For arbitrary time-independent potential flow we actually prove that 
transport is always depleted because of trapping. For time-dependent potential flow or 
generic compressible flow it is shown by explicit examples that transport rates are en- 
hanced or depleted depending on the specific structure of the velocity field. Second, and 
more surprisingly, trapping effects due to the compressibility of the flow can enhance the 
spreading of particles and lead to very efficient ballistic transport. This means that the 
average distance of a particle from its initial position does not vanish, but actually grows 
linearly with time. Note that the small-scale velocity field is supposed to have zero av- 
erage value. Still, because of compressibility, an effective mean velocity emerges in the 
large-scale dynamics. Particles are indeed strongly concentrated or rarefied and sample 
the small scales in a very non-uniform way. It is therefore possible to have an effective 
large-scale mean velocity even if the small-scale flow has zero average value. 

The paper is organized as follows. In Section ^, we formulate the problem and present 
some simple heuristic arguments. The multi-scale technique for ballistic transport is 
presented in Section ^. Simple examples of flows leading to ballistic transport are discussed 
in Section ^. The general conditions ensuring the absence of the ballistic effect, i.e. 
isotropy and/or parity invariance, are finally discussed in Section |^. For the flows where 
the ballistic effect is absent, transport is diffusive in the presence of scale separation. The 
general formalism for the calculation of the effective diffusivities in compressible flow is 
presented in Section ||. The special case of static potential flow is investigated in detail 
in Section |^. The final Section is reserved for conclusions. 
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2 Formulation and heuristics 



Particles advected by a velocity field v and subject to molecular diffusion obey the fol- 
lowing Langevin equation 

^^^ = vix,t) + Wit). (1) 

The function x{a,t) denotes the position at time t of the particle which was initially 
in a. The random process W{t) is Gaussian, independent of v, has zero mean and is 
white-noise in time 

E {Wi{t) Wj{t')) = 2K6ij 6{t - t'). (2) 

The constant k appearing in @) is the molecular diffusivity. The velocity field v is here 
supposed to belong to one of the following classes : (i) deterministic and periodic in space 
and time or periodic in space and time-independent. The period in the various directions 
need not be the same ; (ii) random, homogeneous and stationary or homogeneous and 
time-independent. In both cases the velocity is a prescribed iunciion of a;, and possibly of 
t, and we shall not be concerned with the mechanisms maintaining the fiow. The mean 
value (■) denotes the average over the periodicities for deterministic fiow and the ensemble 
average in the random case. The velocity field has vanishing average value 

{v) = 0. (3) 

Note that this hypothesis does not involve any restriction since we can always perform a 
Galilean transformation to reduce to the case @. 

The Fokker-Planck equation for the probability density Q{x,t) which is associated 
with the Langevin equation (|l|) is (see, e.g., Ref. pO|) : 



dtQ + d-{ve) = Kd^e, (4) 

where d denotes the spatial gradients. The problem discussed in this paper is the large- 
scale behaviour of (^. Specifically, we are interested in the dynamics of the field B on 
scales large compared with those typical of the small-scale dynamics, e.g., the periodicities 
of V. This is the situation encountered in the evolution of spots of particles for large times. 
As it was discussed in the introduction, the goal is to derive effective transport equations 
involving only the large-scale degrees of freedom and to calculate the transport coefficients. 

Before presenting the systematic formalism of the next Section, it is worth to consider 
a simple heuristic argument, in the same spirit as in Ref. for the dynamo problem. 



Taking the average in and trying to write down a closed equation for the mean (B), we 
face the classical closure problem : the Reynolds stress tensor (vB) cannot be expressed 
as a function of (B). The same closure problem appears also analyzing higher-order 
fields correlations p2[, which will not be considered here. However, in the presence of 
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scale separation the situation changes. An approach similar to the Born-Oppenheimer 



approximation in solid-state physics |]23[ becomes feasible. There, the electrons are much 
faster and follow at any moment the ionic configuration. Here, the small-scale dynamics 
has enough time to readjust to the large-scale field and becomes essentially slaved to it. 
It is then possible to express the Reynolds stress tensor by a gradient decomposition as 

J^ = {v,Q) = (H).(e) - Kgv,-(e) + . . . , (5) 

where the remainder includes higher order terms in the gradients. Inserting into (1), 
we obtain 

5t(e) + n-v(e) = «:gv,v,(e). (6) 

This equation describes ballistic transport with mean velocity and diffusion in the 
comoving frame of reference with an effective diffusivity nfy A spot of particles will 
then be rigidly transported by the mean velocity Vb and at the same time its contour is 
deformed, diffusing with diffusion coefficient nfy The corresponding long-time behaviour 
of particles dispersion is 

((x.(t) - x.(0)) {x,{t) - x,(0))) = (H), (H), + 24 1. (7) 

It is remarkable that, despite of the zero average condition (H), a mean velocity Vb appears 
in and the dispersion in (|^ is quadratic in t. The physical mechanisms leading to 
this effect will be made clear by the examples of Section ^. There are nevertheless some 
classes of flow where we can already foresee (the systematic arguments are presented in 
Section ^ that the ballistic term in (^) will be absent and dispersion is linear in t. First 
of all, for incompressible flows the term proportional to the mean (0) should not appear. 
A constant field, being a trivial solution of (^, has in fact no dynamical role. Second, for 
parity-invariant or isotropic flows we also expect Vb = 0. A deterministic parity-invariant 
flow has indeed a center of symmetry, e.g. the origin, so that 

v{-x,t) = -v{x,t). (8) 

For random flow, parity invariance means that the statistical properties of v are left 
invariant by the operations x ^ —x and v ^ —v. Since any vector and its opposite are 
equivalent, it follows Vb = 0. The effective velocity also vanishes for isotropic flow since 
no preferential direction can be picked out. 

Note that two other examples of transport phenomena also leading to a first-order 



dynamics are known : the anisotropic kinetic alpha (AKA) effect and the a-effect 



pT| in magneto-hydrodynamics. For the presence of both these effects the breaking of 
parity-invariance is needed and the AKA effect disappears for isotropic flow as well. The 
main difference with respect to ballistic transport is that the AKA and the a dynamics 
are not purely dispersive and lead to large-scale instabilities in three dimensions. 
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3 Multi-scale theory for ballistic transport 



A general method for dealing with the singular perturbation problems encountered in 
systems with scale separation is provided by multi-scale techniques 0. As discussed in 
the previous Section, we are interested in the dynamics of the large scales. In concrete 
situations the scale separation between small and large scales is obviously always finite. 
The common procedure using asymptotic methods is however to consider the limit of very 
small perturbations, obtain an expansion which is asymptotically valid and use it then 
for finite values of the perturbation (see, e.g., Ref. 0). In our case the units are therefore 
defined in such a way that the typical scales of the velocity field v are 0(1), while the 
large scales where transport takes place are supposed 0(e~^), with e <^ 1. In multi-scale 
techniques a new set of space-time variables X and T (called "slow" ) is introduced. The 
rationale in the choice of the variables is that the large-scale dynamics should take place 
on scales 0(1) in the new variables. For example, since the spatial large scales are 0(e~^), 
we should define X = ex. As about the choice of the time variable, it depends on the 
specific case. Ballistic transport corresponds to a first-order equation in space and time. 
It follows that slow time for this case should be defined as T = et. On the contrary, 
for diffusive transport the rescaling should be T = e^t since the diffusion equation is 
first-order in time, but second-order in space. The different re-scalings illustrate the fact, 
well-known in hydrodynamics, that advection and diffusion take place on different time 
scales. In order to capture both one could use a two-time formalism, as in [^. However, 
for the sake of clarity we prefer treating advection and diffusion separately. In this Section, 
which is devoted to ballistic transport, we define then T = et. The key for overcoming the 
singularity of the perturbation is to pretend that fast and slow variables are independent. 
This allows to correctly capture the dynamics of stirring between large and small-scale 
modes which is missed by regular perturbation theory (see, e.g., Ref. [Q). It follows that 

di^ di + eVi , dt^ dt + edr , (9) 

where we shall denote the derivatives with respect to fast space variables by the symbol 
d and those with respect to slow variables by V. The solution B of (P is then sought as 
a series in e : 

= 0(o)+ee« + e20(2) + ..., (10) 

where all functions depend a priori on both fast and slow variables. Let us then insert the 
expansion ( p!ol) and the derivatives (P) into the original equation (|[). By equating terms 
having equal powers of e, a hierarchy of equations is generated. The general structure of 
the equations is : 

Af^dJ + d-ivf)-Kd'f = g, (11) 

where / is the unknown function and g is known (possibly as a function of the solution 
of lower order equations). Since the operator A has derivatives on the left of all the 
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terms, the solvability conditions {g) = must be satisfied for equation ( pT]) to have a 
solution (Fredholm alternative)^] The solvability conditions fix the dependence of the G*^"^ 
functions in ([10|) on the slow variables averaging over the small-scale degrees of freedom. 
It is therefore not surprising that it is precisely by the solvability conditions that the 
effective equations for the dynamics of the large scales are obtained. 

The equations which are relevant for the analysis of ballistic transport are : 

^e(°) = 0, (12) 

= -9Te(°)-v- ve(°) + 2Ka-ve(°). (13) 

Let us first briefly consider the incompressible case. Since d ■ v = 0, the solution of 
(|12D is a trivial constant field O'-^-* = Q^^\X,T). Inserting this expression into (p!3D and 
using (^ and the periodicity of v (or its homogeneity in the random case), we easily find 
that 

9r(eW) = 0. (14) 

No dynamical process is therefore taking place on time-scales 0(e~^). This is in agreement 
with the known result that (in the presence of scale separation) incompressible flow 
always lead to diffusive transport, which involves time scales 0{e~^). 

The situation changes when a compressible velocity field is considered. Using the 
linearity of (0), its solution 6*^°^ can be expressed as 

e(°) = (e(°))(X,T) m(a;,t), (15) 

where m has unit average value. The solution m is in general a nontrivial function due 
to the interplay between the solenoidal and the potential components of v. By inserting 
(|15D into (0), the following solvability condition is found 

9T(e(°))(x,r) + v(0(°))(x,r) = 0. (i6) 

Equation (|16D is indeed the one describing ballistic transport. The progress with respect 
to the heuristic arguments of Section |^ is that we have now the expression of the effective 
velocity 

Vf,= {mix,t)v{x,t)), (17) 
with the field m being the solution of the auxiliary problem 

dtm + d ■ {v m) — nd"^ m = ; (m) = 1. (18) 

^The solvability condition is, more specifically, that g should be orthogonal to the null space of the 
operator adjoint to A, which is made of constants. Note also that the same multi-scale methods could 
be used to analyze the long-time behaviour of the equation dtO + v ■ d 9 = nd^O. Its adjoint operator is 
closely related to the Fokker-Planck operator (^) . This allows to carry over most multi-scale results valid 
for (0) to the previous equation. 
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Note that (0) has a very simple physical interpretation. The interplay between poten- 
tial and solenoidal components modifies an initially uniform distribution of particles into 
a nontrivial density m. Equation simply illustrates the fact that the small scales are 
sampled non-uniformly with the density m, which is naturally selected by the dynamics. 

In Section ^we shall discuss the conditions ensuring Vb = 0. For generic compressible 
fiow Vb will however not vanish and we have reduced its calculation to the solution of 
(|T8D. Except for a few cases, discussed in the sequel, this equation cannot be solved 
analytically and numerical methods should be used. We shall not dwell here on numerical 
aspects (see, e.g., Refs. |T(]| , [TI|]). We just remark that solving auxiliary problems rather 



than the original equation is very convenient since only small scales are involved in (|T^). 



4 Examples of ballistic transport 
4.1 Flows depending on a single coordinate 

The main features of ballistic transport are well illustrated by the simple class of time- 
independent fiow depending on a single coordinate. For these fiows we show here that the 
effective velocity Vb can be calculated analytically. 
A two-dimensional example is 

v{x)=(^-^{x);w{x)y (19) 

The flow is clearly a superposition of the potential part —d(f)/dx and the solenoidal com- 
ponent w. In the following we shall consider the two-dimensional case for simplicity, but 
our results can be easily extended to more dimensions. The solution of (|T8|) for the flows 
([IqD is the Boltzmann distribution corresponding to the potential : 

rn- r'-rtlrl (20) 

(exp(-0/K)) 

Note that (|^) depends on the potential component only. This property is due to the very 
simple geometry involved in the flows (|T^). It follows from (PO) that the x-component of 



the ballistic velocity vanishes because of homogeneity or periodicity (see eq. (p3D). On 
the other hand, the expression of the y-component is : 

_ {w{x)exp{-(t>{x)/K)) 

which in general does not vanish. Equation (^Tj) provides the analytic expression of the 
effective ballistic velocity. 
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Figure 1: Log-log plot of the dispersions (lower curve) and {y'^if)) corresponding to 

the flow (22), versus time. These curves were obtained by Monte-Carlo simulation. The 
asymptotic slopes are respectively 1 and 2, consistently with the fact that the ballistic 
velocity points in the y-direction. The values for the ballistic velocity generated in the 
simulation agree with the theoretical prediction (23). 



For the sake of concreteness, let us consider in more detail the following simple periodic 
flow : 

V = (sinx; — cosx). (22) 



The averages appearing in (^) can be calculated in terms of the Bessel functions Jq and 
Ii of imaginary argument. The relevant formulae, from are listed for convenience in 
Appendix B. The flnal result is 



(23) 



'y /o(l/«:) + /o(-l/«:) hil/n) 
For large values of the molecular diffusivity k, the Bessel functions can be expanded and 



(H), = ^ + 0(l/«:^ 



(24) 



which agrees with the general expression (|A.7| ) derived in Appendix A. When n is small, 
we can use the asymptotic expansions ( [B.5|) and ( [B.6|) to obtain 



1-- + 0{k') 



(25) 
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Both (^) and ( ^5]) can obviously be also obtained directly from (^Tl) . For small k's the 



field m has a sharp maximum in tt around which it reduces to a Gaussian with variance 
K,. The result (|25|) is then obtained by using Laplace method. The expansion ( p5D is 
in agreement with the fact that ^(vr) = 1 and molecular noise forbids the concentration 
of all particles in the minimum of the potential. A closer look at the streamlines of 
the flow reveals the very simple dynamics associated with the ballistic transport (^). 
The potential component concentrates the particles in the middle of the channel (0,27r). 
Here, the solenoidal part w transports the particles in the positive y-direction. In the 
general case the mechanism is more complicated and cannot be treated analytically but 
this simple example highlights the essential role of both the solenoidal and the potential 
components for ballistic transport. 

4.2 A cellular flow 



Another velocity fleld such that the auxiliary equation (|T8D can also be tackled analytically, 



at least partially, is the following slightly perturbed cellular flow 

V = (cos y + 5 cos x ; sin x) = + x5 cos x. (26) 

The term proportional to the constant 5 is a small compressible perturbation, while for 
5 = the flow (|26| ) reduces to the projection of the incompressible ABC flow (with 
A=0, B=C=1) on the x — y plane. We shall derive here the asymptotic expression of the 
effective velocity valid for large Peclet numbers. 

The incompressible BC flow consists of square convective cells with separatrices 
parallel to the diagonals. The mechanism of effective diffusion at high Peclet numbers 
for this type of flow is well understood p, P). Particles within the cells simply circulate 
along the streamlines and their density rapidly becomes almost uniform. On the contrary, 
particles close to the separatrices may cross the cell boundary because of molecular noise. 
The width of the layers where these transitions typically occur is proportional to y/n. 
For small /t, it is therefore possible to calculate the asymptotic behaviour of the effective 
diffusivity by using boundary layer theory. The stream function of our BC flow 
coincides with the one considered in Ref. |]^ with the transformations 

/ x + y , x-y TT 

and u = ^/2, d = V2n, P = 1 (see eqs. (2) and (3) in Ref. 0). The corresponding 
expression of the effective diffusivity derived in Ref. ^ (eqs. (26) and (28)) is: 

K^" '^k'/\ (28) 
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Figure 2: Contour lines of the function m corresponding to the shghtly perturbed flow 
(26), obtained by solving the auxiliary problem (18) numerically. The parameter values 
are S = 0.05 and k, = 0.01. 



On the other hand, the effective diffusivity for the incompressible BC flow can be 
expressed as (see, e.g., Ref. [|lOl) : 



K,6, 



1 



Q/3 



where s has zero average value and satisfies 

dts+ (v^ ■ d) s — kO^ s = —v^. 



(29) 



(30) 



The effective diffusivity, which is in general a second-order tensor, reduces for the BC flow 
to a scalar thanks to the following symmetries : 



X y 



TT 



y i-^ X 



TT 



(31) 



and 



xi-^vr — x; y ^ y — Ti \ vl.^ —v^ 



X ) 



(32) 



An immediate consequence of (|3T|) is that the diagonal components of the tensor in ( pOf ) 
are equal and the symmetry (^) implies that the non-diagonal components vanish. It 
follows that 

< = -K^.), (33) 
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Figure 3: Same as in Figure 2, with parameters 6 = 0.05 and k = 0.001. 



where Sy is the zero-average solution of the equation 



dt Sy + (^v^ ■ Sy — K d'^ Sy = — slu 



X. 



(34) 



The crucial remark is now that equation (|3^) coincides with the auxiliary equation 



([TsD for the flow (|26|) at the dominant order in 5. When the solution of (pTSD is sought as 

m = m(°)+5m(^)+52m(2) + ..., (35) 



it is easy to check that m*^'') = 1 and m^^^ = —Sy. It follows from (17), (2S) and (0) that 
the asymptotic behaviour of (Vfe)y is 



25 



O.7490F 



(36) 



Note that the behaviour in k of terms of higher order in 5 is not under our control in this 
expansion. 

The x-component of the ballistic velocity vanishes for arbitrary 5 and k. Independently 
of the value of 5, the flow (EB) has in fact the symmetry (R2|) : 



By using (|37|) in (0), we immediately obtain 

m{7i — X ; y — n) = m{x ; y), 



Vy ^ Vy. 



(37) 



(3^ 
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and 



(HL = (t;.m) = 0. 



(39) 



The prediction ( p6D is compared in Table 1 to the values obtained by solving numeri- 
cally the auxiliary equation (|18D for the flow (^). The equation has been solved by using 
pseudo-spectral methods ||2^ and the parameter 6 = 0.05. The contour plots of the field 
m are shown in Figs. |^ and |^ for k = 0.01 and k, = 0.001, respectively. The resolution is 
256 X 256, which ensures that the boundary layers are properly resolved. 





Spectral 642 


Asympt. Pred. 








1. X 10-1 
5. X 10-2 

1. X 10-2 

5. X 10-3 


2.0 X 10-2 
1.5 X 10-2 

7.1 X 10-3 

5.2 X 10-3 


2.4 X 10-2 
1.7 X 10-2 

7.5 X 10-3 
5.3 X 10-3 



Table 1: The ?/-component of the ballistic velocity (H)^ for the flow (^) with S = 0.05. 

It is evident that the quality of the prediction (^) improves when k is decreased. This 
is in agreement with the asymptotic nature of the boundary- layer technique which has 
been used in the calculations. 



4.3 Random flows with short correlation times 

It is known that the effective diffusivity for incompressible flow with short-correlation 
times can be calculated exactly (see, e.g., Ref. fl^])- The aim of this Section is to show 
that short correlation times also allow to calculate the effective ballistic velocity for com- 
pressible flow. 

Specifically, the velocity field is supposed to be a homogeneous, Gaussian random 
process having zero average value and correlation function 

Mx, t) vp{0, t')) = F^p{x) 5{t - t'). (40) 

In order to calculate the average (|T^ we can now use the formula of Gaussian inte- 
gration by parts (see, e.g., Refs. p^, p8[). The average is then expressed as 



V, = {mv) = J dx'dt'{v{x,t)v{x',t')){^^0^^). (41) 
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From the auxiliary eq. ([I8| ) and the (5-correlation in time it immediately follows 



V, = - {v{x,t)id-v){x,0))dt, (42) 
Jo 

where the stochastic differential equation (ITS) is interpreted a la Stratonovich since we 



are interested in the physical limit of short correlation times. Eq. ( ^2]) had already been 
derived in Ref. |^ using path-integral methods. Note that the average appearing in ( ^^ 
vanishes for incompressible flow. It also vanishes for parity-invariant and isotropic flow, 
in agreement with the results of the next Section. 

An alternative derivation of (^) is to neglect the diffusive term in the auxiliary equa- 
tion ([18|) and integrate on the characteristics : 

exp(-f_^{d-v){x{a,s);s)ds) 

m{x,t) = ) (43) 

{exp [- J_^{d ■ v){x{a, s) ] s) ds)) 

The characteristics are the Lagrangian trajectories defined by the equation 

^^^ = v{x{a,t)-t), (44) 

where a is the initial position of the particle. In (^31) the integration has to made on 
the trajectory such that x{a,t) = x. Inserting ( ^5D into (^) we obtain an expression of 
the ballistic velocity as time integral of Lagrangian correlations. This is the equivalent 
of Taylor's formula for effective diffusivities. The presence of the exponential makes 
it however very difficult to use. When the flow has a short correlation time one can 
however simplify ( ^3D by exploiting the fact that Lagrangian and Eulerian statistics tend 
to coincide. Using again Gaussian integration by parts (^2]) immediately follows. 



5 General conditions for the absence of ballistic 
transport 

The aim of this Section is to identify some general classes of flows where ballistic transport 
is absent and pure diffusion takes place. 

A first class where Vb = is the one of incompressible flow. This has already been 
shown by (|1^ in Section ^ 

Let us then consider parity-invariant flow, i.e. velocity flelds having at least one 
center of symmetry according to the deflnition (|^) in Section ^. It follows from (^ that 



the solution m of (|T8D having unit average value satisfles 

m{—x,t)=m{x,t). (45) 
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Figure 4: The dispersions (lower curve) and (y'^it)) for the flow (47). The straight 

line is proportional to t. 



In the random case, the statistics is invariant under the operations 



X 



-x: 



V 



-v; 



m 



m. 



(46) 



An immediate consequence of (|^) and (45) (or (46) for the random case) is the vanishing 
of the ballistic velocity Vb = 0. 

To give a direct confirmation of these considerations we present in Fig. ^the dispersions 
for the flow 

(sinx ; sinx + sinSx). (47) 



V 



This flow depends on a single coordinate as discussed in Section |4.1| . Unlike (|22D , the flow 
(1^ has however the center of symmetry x = vr. It follows then that ballistic transport 
is forbidden. The difference is best appreciated by comparing Figs. 1 and |[ In the flrst 
case {y'^it)) oc while in Fig. ^ {y'^if)) oc t, corresponding to diffusive transport. 

We consider now homogeneous and isotropic random flow. These flows are deflned, 
According to the deflnition in pDI], these flows are such that their statistical properties 



are unaffected by translations and/or rotations accompanied by a simultaneous rotation 
of the laboratory frame. Note that, unlike the deflnition in Ref. |29|, we distinguish here 
between isotropy and statistical parity-invariance, e.g. a 3D helical flow can be isotropic 
even if it is not parity- invariant. Since no preferential direction but v appears in (|l^), 
the statistical properties of the couple (m, v) are invariant under rotations. Homogeneity 
and the absence of a preferential direction imply then that the single-point correlation 



Vb = {v{x,t)m{x,t)) = 0. 



(48) 
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In practice, information on isotropy or parity-invariance of the flow may not be avail- 
able. To check whether baUistic transport is present one could then use the expansion in 
powers of ^ derived in Appendix B. To leading order, this expansion shows that 

V, = --■ {vH-' {d ■v))+ O , (49) 



where Ti. ^ = {dr — d"^) ^ is the Green's function of the heat diffusion operator. For 
time-independent flows, the latter expression reduces to 

V, = --■{v'<P)+ O(^) , (50) 



where represents the incompressible component of the velocity field and the potential 
of the irrotational component. The averages in the asymptotic expressions (^) or ( ^0|) 
are easy to measure and can be used to estimate whether the ballistic drift is present or 
not. 

Finally, we investigate time-independent potential velocity fields, i.e. satisfying 

v{x) = —d(j){x), (51) 
where is the potential. The stationary solution of (|18]) for this class of flows is the 



Boltzmann distribution (see, e.g., Ref. pO[]) 



It follows from periodicity (or homogeneity) that 



exp f-^) 

^ = r^- (52) 



Vb oc (i;exp -- ) = exp -- ) =0. (53) 



K 




Ballistic effect is therefore impossible for time-independent potential flow. Note that the 
condition of time-independence is essential for this result. Indeed, if depends on time 



the solution of (|I8D is not the Boltzmann distribution. Furthermore, let us consider the 
one-dimensional flow 

v{x,t) = cos{x — Kut) + acos{x + nut). (54) 

Using (|A.5| ) in Appendix A, it is easy to check that the first coefficient of the perturbative 
expansion in 1/k of the ballistic velocity is 

V^-(r^'-^)^AOa)^ (55) 



The fact that (|55|) does not vanish illustrates the fact that ballistic transport is possible 
for time-dependent potential flow. One can also easily check that the expression (^), 
valid for compressible (5-correlated flow, does not generally vanishes for potential flow. 



16 



6 Effective difFusivities for compressible flow 



In the previous Section we have considered advective effects in large-scale transport by 
compressible flow. It is physically clear that diffusive effects will also be present. In 
order to systematically treat them, multi-scale techniques can again be used. To capture 
both advective and diffusive effects, two independent slow times Ti = et and T2 = e^t 
are needed. The solvability condition at order e gives the advective effects, as shown in 
Section ^, while at order we also include diffusive effects (see Ref. P). For the flows 
where the ballistic velocity Vb vanishes, no advective time is needed. The formalism is 
then the same as in Section ^, with the only difference that the slow time is now defined 
as T = eH. Since there are just technical differences between the two cases Vb 7^ and 
Vb = 0, we shall treat in detail the latter and give only the final results for the former. 

As in Section |^ we denote the fast variables by x, t and the slow space variable is 
defined as X = ex. The analogous of equation (P) for the derivatives is now : 

d,^di + eV, , dt^dt + t^dT , (56) 

The solution G of the Fokker-Planck equation (^ is again sought in the form of a series 
in e. When (^6]) and (|l^) are inserted into (^), a hierarchy of equations of the form (|lT]) 
is recovered. For the analysis of turbulent diffusion the following three equations are 
needed : 

^e(°) = 0, (57) 

^40(1) = -t; . ve(°) + • ve(°) , (58) 
^0(2) _ _5^e(o)_^. v@{i) + 2/€5- ve(') + fi:V'e(°). (59) 

The operator A is defined by (0). Using the linearity of (0), we can express the solution 



as in (p!5|), where the field m satisfies (|T8|). We can then insert ([15|) into (pS]). Using again 



the linearity of the operator A, we obtain: 

e« = x- v(e(°)) + m(eW). (60) 

where the vector field x satisfies 

dtX + 9 ■ {vx) ~ ^ = m + 2k d m. (61) 

The underline in the l.h.s. of (^TJ) is meant to make it clear that in the scalar product the 
components of d have to be contracted only with those of v, and not of X- As expected, 
the solvability condition for (|6TD requires Vb = (which we have supposed to be true). 
In the cases when Vb does not vanish, the calculations with the two-times formalism give 
the following modified equation 

dtX + 9 ■ (vx) — K, d'^x = ('^ ~ ^b) + 2Kdm, (62) 
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where Vh is given by (|T^. Eq. (|6T|) is clearly a particular case of (|62D . 

The equation for the large scales dynamics emerges as the solvability condition for 

m 

ar(eW) = «:f^v.v^(e(°)), (63) 

where the effective diffusivity tensor is 

f^al3 = l^^c^P - \ {{VaXn) + • (64) 

When Vb 7^ 0, the effective equation for the mean value of the field becomes 

Q^iem) + V, ■ v(eW) = K%s/^s/p{Q'^'^), (65) 



where the effective diffusivity has the expression ( |6^ and the x field is the solution of 

The diffusivity tensor /tf^ might of course be anisotropic, reflecting the fact that the 
flow could be non-invariant under arbitrary rotations. The term in the square parentheses 
in (|6^ does not have a definite sign for generic compressible flow. The quantity — k 
is in fact known to be positive for incompressible flow and in the next Section we show 
that it is always negative for static potential flow (see also the following example). 

There are some special instances where the effective diffusivity can be calculated ana- 
lytically, e.g. the flows of the class (p^ . We shall not dwell on the details here since the 
solution of (|TB|) is the Boltzmann distribution ( pOD and the solution of (^Tj) is essentially 
the same as the one which will be discussed in Section [7.2| . The final result is that is 
given by (BDI), the non-diagonal components vanish and 



2 



which is clearly larger than k. The function is defined by 



■\\){x) = J w{z) exp y —j dz, (67) 

and it is still a periodic function since the ballistic velocity is supposed to vanish. 



For generic compressible flow, the solutions of the auxiliary equations (18) and 



are not generally known and one should then resort to numerical methods, as in the 



incompressible case [10, 11 
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7 Static potential flow 



One of the results derived in Section ^ is that no balhstic transport is possible for static 
potential flow. In the presence of scale separation, transport by these flows is therefore 
always diffusive. This is rather intuitive from a physical point of view. Particles are indeed 
concentrated into the minima of the potential where the velocity however vanishes. For 



the flow in Section [4.1| one immediately realizes, for example, that without the streaming 
in the ^-direction due to the solenoidal part no ballistic advection would be possible. As 
a matter of fact, we expect that even diffusion should be very sensitive to the presence 
of potential wells. Let us indeed consider a stagnation point a^o {v{xo) = 0) and the 
quantity 

r = v{x) ■ (x - xo). (68) 
When T is expanded in a Taylor series, the first non-zero term is generically quadratic : 

r ^ S^,ix - xo). {X - xo), = \ {d.v, + d,v.) {X - xo). {X - X,), (69) 

In order to have strong trapping the velocity field should be such that the particles are 
driven back to a^o; i-e- ^ < everywhere in a neighbourhood of xq. This is by definition 
what happens in the minima of the potential for the flows (|5T|). On the contrary, this 
strong form of trapping is impossible for incompressible flow. The tensor S is in fact 
traceless and not all its eigenvalues can be negative[[ The conclusion of these simple 
heuristic arguments is that a depletion of transport is expected for static potential flow. 
This expectation is confirmed by the systematic results of the next Section showing that 
the effective diffusivity is indeed smaller than the molecular diffusivity. 

7.1 Depletion of transport 

For static potential velocity fields of the class (|5l|) the calculation of the effective diffusivity 
can be reduced to the solution of one auxiliary equation only. We can then show that 
diffusive transport is always depleted. 

The solution of the first auxiliary equation (|57| ) is the Boltzmann distribution : 



e(o) ^ (e{o))(x,r) """"^^ . (70) 



Inserting (170) into (|6TD, we obtain: 



dt'E + d-{v'E) - Kd'^ = exp -- , (71) 




^Note that a weak form of trapping is however possible 
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where S = x(exp (~^)) and has zero average value. It is then convenient to define 

3. (72) 




The equation for T is easily derived from ( [7TD 

dtr-{v-d)r-Kd^r = v. (73) 

We can now recast the averages appearing in the expression of the effective diffusivity 
in a more convenient form. Let us indeed multiply the components a and (3 of (|73|) 
by Ejs and Sq, respectively. By summing, taking the average and using the definitions of 
H and T, we obtain after some algebra: 

^ (vaXp) + (v^Xa) = -«^(exp ( ) {dT^) ■ {dT,)) /(exp [~^]). (74) 



2 



K \ K 



In order to derive this equality we have used the periodicity (or homogeneity) of f . It is 
evident from ([7^) that the correction to the molecular value which appears on the r.h.s. 
of (|^) is negative definite. This proves that transport is always depleted. 

Let us then show that, despite of the fact that the correction due to the fiow v is 
negative, the whole effective diffusivity /tf^ is still positive, i.e. the absolute value of the 
correction is always smaller than k. Inserting the expression of x iii terms of T into 
and integrating by parts we have 



(exp(-fj) 

where n is a generic direction and d\\ = {n ■ d). Using the Schwartz inequality leads to 



(exp {-i) 


(a||(T-n))2)V2 


(exp (-!) 


)l/2 



(exp (-^) 


(9(T-n))2)l/2 


(exp 




)l/2 



Equation ( |7iD has been used to derive the last equality. The consequence of (^) is that 
the effective diffusivity tensor K^jg is indeed positive definite, in agreement with the known 
fact that the original Fokker-Planck equation is stable . 

We have then shown that a static potential flow always leads to depletion of transport. 
This should be contrasted to the case of incompressible flow, where transport properties 
are always enhanced. Remark that for time- dependent potential flow the result just proved 
does not generally hold. The time-dependence can in fact destroy the trapping effects due 



20 



to potential wells and lead to an enhancement of transport. One can easily provide an 
explicit example by taking again the flow with a = 1 (the flow is then parity-invariant 
and no ballistic effect is present). The analytic expression of the effective diffusivity is 
not known. For our aims it is however enough to perform an expansion for large k quite 
similar to the one presented in Appendix A. The first term of the expansion of is 
obtained by some simple algebra : 

n^ = n + -^4i^2+Oil/n'). (77) 

K (C(J^ + 1)^ 

The fact that can be either larger or smaller than k indicates that no general rule 
can be expected. Enhancement or depletion of transport properties are both possible, 
depending on the detailed structure of the flow. 



7.2 The one-dimensional case 

The calculation of the effective diffusivity has been reduced to the solution of the auxiliary 
equation ([7T|). While in the multi-dimensional case no general solution is known, the 
equation can be solved analytically in ID. The reason for this simplification is that in 
one dimension the order of the equation can be lowered, thus reducing it to a first-order 
differential equation. We shall consider in detail the case of periodic, deterministic flow. 
The random case is handled similarly. 
It follows from equation ([73|) that 

A(T + .) = .expQ, (78) 

where c is a constant and x is the spatial coordinate. One more integration leads to 

T = c J exp ^— ^ dy + c' — x. (79) 

The constants c and c' are fixed by imposing the conditions that the function x appearing 
in ( [72|) is periodic and has zero average. For the calculation of k,^ only the value of c 
turns out to be actually needed. By imposing the periodicity of (|79|), we finally obtain 
the expression of the effective diffusivity : 

^ ^ (80) 

(exp(^))(exp(-^)) ^ ^ 

It follows from (PI1| ) that for small values of the molecular diffusivity the ratio k^/k 
is exponentially small and has the structure of an Arrhenius factor. The term in ( pOD 
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is indeed proportional to the probabihty for a particle to jump out of a well of depth 
equal to the difference between the absolute maximum and minimum of the potential 
(see e.g. Section 5.10 in Ref. |2^). The mechanism of transport is thus essentially a 
random walk between the minima of the potential. Particles are trapped for very long 
times in the bottom of potential wells. Because of large fluctuations in the noise, they can 
occasionally jump out of a minimum falling into the adjacent one. The same mechanism 
evidently works also in the multidimensional case and is responsible for the depletion of 



transport found in Section 7.1 



For Gaussian random flow the averages in (^) can be easily calculated : 



exp 



Note that for transport to be diffusive at sufficiently large scales, the condition 

(0') < OO, (82) 



must be satisfied. When ( |52D is not satisfied, the hypothesis of scale separation breaks 
down. The small-scale typical time for relaxation to local equilibrium is indeed of the order 
of the typical Arrhenius time needed to jump out of the wells. When (|^) is violated, 
this time is divergent because wells of arbitrary depth are present. There is therefore no 
time which is much larger than the typical times of the small-scale dynamics. Standard 
diffusion is never observed and transport is sub-diffusive. An explicit example is provided 
by the case when is a Brownian motion. Ya. Sinai found in Refs. 0, ^ that in this 
case the dispersion does not vary linearly with time but 

((x(a,t) -x(a,0))2) ~ (logt)^ (83) 



Condition (|82|) is not satisfied also when is a fractional Brownian motion. This problem 



together with other examples of anomalous subdiffusive behavior and applications to 



statistical mechanics systems are discussed in detail in Ref. p3 . 



8 Conclusions 

Turbulent transport by compressible flow has been investigated. For generic compress- 
ible flow having non-vanishing solenoidal and potential components, ballistic transport 
typically takes place. For isotropic and/or parity-invariant flow the ballistic velocity van- 
ishes and transport is diffusive. The properties of the effective diffusivities are strongly 
dependent on the compressibility of the flow. For incompressible flow it is known that 
diffusion is always enhanced in the presence of small-scale velocity flelds. On the contrary. 
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we show that for static potential flow the effective diffusivities are reduced with respect 
to their molecular value. Trapping effects due to potential wells are responsible for this 
depletion. For time-dependent potential flow and generic compressible flow, depletion or 
enhancement of transport are both possible. No general rule can be drawn and the type 
of transport depends on the precise structure of the flow. An essential role is in particular 
played by the geometry of channels and stagnation points of the velocity field. 
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Appendices 



A Expansion of the ballistic velocity for large 
molecular diffusivities 

In the limit of large molecular diffusivities k, the ballistic velocity can be expressed as a 
perturbative series 

oo (n) 
n=l ^ 

The vector coefficients c^""^ are defined by 

cW = (t;m(")), (A.2) 
where the functions m*^"^ are obtained recursively as 

^(n+l) _ I^^^(n)^ n>l, (A.3) 

and m^^^ = 1. The operator Ti, is the heat operator 

']-i = dr-d\ (A.4) 

and r = /tt. The recursion relation (|A.3|) is simply obtained by applying the inverse of the 
heat operator to (|THp and expanding in powers of Note that for periodic, deterministic 
flow no problem in the inversion of Ti is encountered. The kernel of the operator is in fact 
made up of constants which however cannot show up in the r.h.s. of ( [A. 3D because of the 
divergence operator. In the random case the convergence of the integrals involved in the 
coefficients c*^"-* depends on the specific behaviour of the correlation functions of v. The 
first coefficient in can be easily calculated 

c^^^ = -{vH-^d-v)). (A.5) 

The velocity field v can be generally decomposed as the sum of a potential and a solenoidal 
part 

v = -d(l) + v^, (A.6) 

where is divergence-free. For time-independent fiow, (|A.5|) can be then reduced to the 
interesting form 

(A.7) 



c 



(1) 



which illustrates the importance of correlations between the potential and the solenoidal 
components of the velocity field for ballistic transport. 



24 



B Some useful formulae 



In this Appendix we list the formulae which are needed for the calculation in Section |^ of 
the ballistic velocity for the flow (p2]). The formulae are taken from 



(B.l) 



IJx) ~ 



exp(— /3cosx) cosxdx = m Ji{i[3) ; 

exp(— /5 cosx) dx = Ti Jo{iP) ; 
Jn{iz) = i^In{z) ; 

{-If r(n + A; + l/2) 



E 



fct^o (23^)"^ ^' T{n-k + 1/2) 



-X] 



V2 



E 

fc=0 



(-l)^' r(n + A; + l/2) 
(2x)'= k\V{n-k + 1/2) 



(B.2) 
(B.3) 

(B.4) 
(B.5) 
(B.6) 



The functions J„ and /„ denote the Bessel functions of real and imaginary arguments, 
respectively. In both ( [B.4| ) and ( [B.5| ) the variable x is supposed to be real, large and 
positive. The two couples of formulae ( p.l|) , ( p.2|) and ( |B.5|) , ( p.6| ) correspond to 3.915 
and 8.451 in [H. The O and (piP are 8.406 and 8.445. 
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